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Abstract: 


In highly stratified estuarine waters, with significant sediment loads, it is fairly difficult 
to predict the 7,,, the time scale that characterizes the bacterial inactivation or death, 


using mathematical formulae, so that suitable experiments must be done using tracers to 
measure physical dilution and microbiological methods to measure bacterial 
disappearance or inactivation. The location and characteristics of a second sewage outfall 
for the city of Montevideo is being determined. In this framework, several experiments 
with tracers were done to measure the self-purifying capacity of the coastal waters at the 
selected site. Due to the difficulties in the design and interpretation of the measurements, 
special mathematical models were constructed and applied to help in obtaining a deeper 
understanding of the complex advection-dispersion-inactivation processes. In this report 
the construction and the solution of the equations of an analytical model that guided all the 
subsequent modelling activities is described. The tracer and the bacteria fields are 
obtained as products of a function of time and the horizontal coordinates by a series of 
decaying exponential functions of time weighted by functions of depth. The quotient 
between the bacteria concentration and the tracer concentration, that is commonly used to 
determine a global value of the 7,,, taking suitable average values after a long enough 
transient, is studied. A procedure that may be used to summarize the field of local 
inactivation times by a few parameters, determined by tracer experiments and 
microbiological measurements, is discussed. A half-empirical correlation formula 
between the 7;, parameter and the main variables of the body of receiving waters is 
suggested from the results of tracer experiments and other sources of information. The 
advantages and limitations of the modelling approach used in this case are assessed and 
some possible improvements are briefly outlined. 
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(A) Introduction: 


When sewage is discharged in the coastal waters of a river, an estuary or the sea, the 
bacterial concentration decreases due both to physical dilution and to death by adverse 
conditions in the body of receiving waters. The killing rate of micro-organisms due to 
the body of receiving waters depends of temperature, salinity and ultraviolet radiation, 
amongst other environmental parameters, such as other pollutants, nutrients and possible 
predators. 

The inactivation or death is expressed as a 7) value, which is the time it takes for 


bacterial concentration to decrease by a factor of 10 beyond the decrease due to dilution. 
As the nature of this process of self-purification is not completely known, it is in general 
not possible to predict the value of 7, from the properties of the receiving water, albeit 


there are some half-empirical correlations that may be used to obtain preliminary 
estimations. So, an experimental determination of 7,, must be done, as a rule. 


To do this, a population of micro-organisms is injected in the site and its concentration 
is measured at different instants and in different positions in the body of receiving 
waters. 

To distinguish the effect due to advection and dispersion of the population, from the 
killing effects of the environmental factors, a tracer is usually used. The tracer is 
injected jointly with the population of micro-organisms, and its concentrations are 
measured at the same instants and positions. Then, if the tracer behaves as a 
conservative pollutant, and if it is transported in the same way as the bacterial 
population, its dilution must be due only to transport processes. If we use a radioactive 
tracer with a half-life short enough to be taken into account, it is always possible to make 
a correction for tracer disintegration and after doing that, continue as it were a 
conservative tracer. In this case the time evolution of the quotient between the 
concentration of still alive bacteria and the tracer concentration (corrected for 
disintegration if necessary) measured at the same instants and at the same positions, 
would give us a measure of the capacity for self-purification of the receiving waters. 
With this information, together with information about currents and other data, that are 
needed as inputs to be introduced in suitable mathematical models, it is possible to 
calculate, as outputs of these models, what would be the micro-organism concentrations 
in certain critical points (near shorelines, for example) given a location, length, discharge 
flow, and so on, of an intended outfall discharge of sewage waters. 

In estuaries like Rio de la Plata, it is fairly common, mainly during summer, early 
autumn and late spring time, to find a strong vertical stratification. The salinity (and 
density) increases significantly from the free surface to the bottom, and often two or 
three layers of different densities and salinities can be observed. 

Temperature differences between surface and bottom exist, but are usually less 
significant. 

A variable charge of suspended sediments blocks the ultraviolet rays, whose killing 
effect is thus produced mainly in a layer of water adjacent to the free surface, and much 
less or not at all in grater depths, near the bottom. 


As a consequence of all these variations, the kinetics of bacterial death is a fairly 
complex and depth dependent process. So, it is reasonable to ask about the meaning of 
the experimental results obtained during the above mentioned experiments intended to 
quantify the self-purifying capacity of estuarine waters. We will discuss this meaning 
using a fairly simple mathematical model of the processes of advection, dispersion and 
inactivation. 


(B) A simple transport model with decay terms dependent of depth: 


We start from the following field equation for the concentration c(t, x, y, z)of a non 
conservative pollutant: - =L,|c]+L,|c]- K(z)c [1] 
The operator L, describes the horizontal advection and turbulent diffusion, while the 


operator L, describes the vertical turbulent diffusion, and if it is significant, the 


sedimentation process. Both operators, applied to the concentration field, are linear. 
The hydrodynamics related with the advection processes is assumed to be given. 
At this level of description the horizontal eddy diffusivities may be scale dependent. 
However we suppose that the field of local inactivation rates K (z) depends only of the 
vertical coordinate z . (The z axis points upwards from the bottom towards the surface of 
free water). This dependence allows us to take into account the vertical salinity gradients, 
the varying effects of ultraviolet light due to the absorption processes, and so on. 
If E,(z) is the vertical eddy diffusivity and w, is the sedimentation velocity, we can pose 
the following boundary conditions (Nihoul, 1975, chapter1): 

Oc 


At the surface: - E,, wc =0 [2a]. Atthe bottom: -E£, ae wae viel “| [2b] 
iz ve 


In [2b], zc, is the shear stress at the bottom, while 7, is a critical value of zc, such that 
above it the erosion rate is greater than the deposition rate. This is only an approximate 
boundary condition (Fredsoe and Deigaard, 1992) but it suffices for our purpose. 

Now let us consider a layer of the receiving water body of depth H (in our case less than 
10 meters) and several hundred meters of horizontal extension/ . Even if the vertical eddy 
diffusion coefficient E, is one or two orders of magnitude less than the horizontal eddy 


diffusivities E,, the vertical time scale of diffusion +, = Lyf is at least an order of 


magnitude less than the horizontal time scale of diffusion t, = "Ye and the time scale of 
h 


sedimentation tr, = iE We may further assume that 1, is at least an order of magnitude 
Ww 


less than the time scales of the tidal currents and the diurnal variations in the light 
intensity that reaches the surface of the water body. 

Although the horizontal operator may show an explicit dependence of time and 
horizontal coordinates, (due to the advection terms and the scale effect on the horizontal 
eddy diffusivities), the operator L, and the inactivation constant K may be considered as 
independent of time when trying to solve Equation [1]. 

Furthermore, to construct a simple mathematical model, we suppose that L, is 
independent of the horizontal coordinates. 


Then, substituting the ansatz c(t,x,y,z)=y(t,x,y)Z(z) in this equation and introducing a 
separation of variables constant k? we obtain: 


~i,[z]+ K(2)z =k2.Z [3] = Sa yew [4] 


Introducing the auxiliary field y9(v, x, y)=e* ‘v(t.x, y), it verifies the horizontal advection- 
diffusion equation like an already vertically well mixed and conservative pollutant or 


kt 


tracer: “Ps = L, [vo] Meanwhile y is given by: wewee 
The vertical operator is given by L, [d--2[e, ()Zc}-w, “ so equation [3] may be 
Iz, Iz, iz 
written: 
4 [206 Lale)|-m, Beds xtaizte)=4220 5 
dz dz dz 


Ifw, =0 , the differential operator is already in self-adjoint form. 
In any case, it can be put into a self-adjoint form multiplying the equation by the 
. < dz 
function: Plz ~exf + Ww, al [6] 
rere wel) 
From the boundary conditions [2a], [2b] applied to the ansatz, we have, if 
z=Orepresents the bottom and z=H the surface of the water body: 


£, (0) »,2(0)- v. 2001-2 [7a] (1) Ay Z(H)=0 [7b] 


Ss 
he 


From equation [5] in self-adjoint form, jointly with the homogeneous boundary 
conditions [7], we have an eigenfunction (Z(z))-eigenvalue ( k* ) problem of the Sturm- 
Liouville type. If £,(z) is positive and K(z) non-negative, and both functions are 
regular enough, there is a denumerable set of eigenvalues k,*: 0 <ky’ < ky < ky’ <...< 


k,°<... There is a single eigenfunction Z,(z) corresponding to each eigenvalue k,’, and 
the eigenfunctions corresponding to different eigenvalues are orthogonal relative to the 


weigth p(z): eG), (z)Z,(z)dz =0 Pp¥q 


H 
(From now on we consider normalized eigenfunctions: | p(z\Z,(z))7dz=1). 
0 


The initial condition for the concentration field depends of how the injection of the 
pollutant or tracer is made. 

Usually it is possible to pose: c(O; x, y, z)=M,.6(x)d(y) f(z) [8] 
Here M,is the total mass of tracer or pollutant suddenly injected, while 5(x)and 6(y)are 
Dirac’s distributions along the horizontal directions. They stem from the assumption that 
the horizontal size of the mass of tracer or pollutant, immediately after the injection, is 
negligible in comparison with the horizontal scales of transport by advection and 
turbulent diffusion. 

The function f(z)reflects the vertical distribution of the injected tracer or pollutant 


H 
immediately after the injection. It verifies [ f(z)dz=1. Here z, gives the vertical 
: 


coordinate of the deepest point of the injection, so H-z; is the depth of injection 
measured from the surface. 


However, for an asymptotic analysis of the process, it is possible to neglect the details of 
the vertical distribution of the injection. 
The solution of [1] with boundary conditions [2] and initial condition [8] can be 


represented as: c(t,x,y,z)= Letty, (tx, y)Z,,(z) [9] 
n=0 
Then: c(0, x, y,z)= ¥ y,(0;x,y}Z,,(z) 
n=0 
AH 
From [7] it follows: y,(0;x, y)=Mo {fate (2)Z, le alsa) [10] 
0 


This is the initial condition that we must use jointly with [4] to findy, (tx, y). 

These functions are solutions of the same advection-diffusive transport equation 
(Equation [4]) jointly with the same boundary conditions for each of them. 

The only difference appears in the initial condition for eachn , given by [10]. 

So, if yo(t;x, y)is the solution of [4] with the same boundary conditions (not yet specified) 
but with the initial condition y,(0;x, y)= 6(x)6(y), we have: 


v,(t:x%,y)=Mof,Woltixy) [11a] tn SO OVe (zldz [1 1b] 
0 


The only difference between the different functions y, appears in a proportionality 

factor, but the space-time dependence, as given by y,(t;x, y), should be the same for all of 

them. As a consequence, [8] may be written thus: — c(t; x, y,z)=Mp.g(t.z)yo(tsx.y) [12] 

Here, by definition, alt.z)= Dee" f, Z,(2) [13] 
n=0 

From [12] we can calculate the concentration field of the pollutant (or tracer) for t>0 and 

for each point of the body of estuarine waters that we are studying. 

So, assuming that we can determine k,,” and Z,,(z), the only task that remains is to give a 
concrete solution y(t;x,y) for the horizontal transport. But as we are going to consider 
the quotient between the concentration field c,(t,x,y,z)of a non-conservative 
pollutant and the concentration field for the tracer (considered as 
conservative) c,-(t;x, y,z), it is not necessary to derive a formula for y, (t; a y). 

The reason is that y,,(¢;.x, y)must be the same for both, pollutant and tracer. 

The components of the advection field are the same, the turbulent diffusion coefficients 
are the same, the boundary conditions fory, (t; bcs y)are the same and the initial 
conditions fory, (t; x, y) are the same, so the field y, (es x; y) must be the same for both 
of them. 


(C) Calculation of eigenvalues and eigenfunctions by a regular perturbation 
approach 


Let us consider now the eigenvalues k,’ and their corresponding eigenfunctions Z,, (z) ‘ 
From equations [5], [6] and [7], following a standard procedure we derive the relation: 


— jac (e{ 02) +24) LW, ( £.7,2(0) alnz,e) [14] 


dz c 


It is possible to obtain fairly good approximations to the kinetic coefficients substituting 
the exact eigenfunctions by approximate ones, derived by regular perturbation theory. 
The vertical sedimentation-diffusion equation can be written in self-adjoint form: 


-ilzleoxkzZ=e@ pz [15] , é£z]= 4 p62 | with p given by [6]. 
x 


H H 
Let us introduce the averages over depth: Ky = — [K(z)dz [16] Ey = — fé,(z)az [17] 
0 0 
If we substitute, in equation [5], K(z) and E,(z) by their average values, it is fairly easy to 
obtain a closed form solution for the eigenfunctions and the eigenvalues. 


When diffusion time +,=°/ is less than one third of the sedimentation timer, = Ea , 
Y EE i 


Ss 


the eigenfunctions are given by: 


ty 
Zonl2)= | 2Cye ts " cof “12, n=0,1,2,3.... [18] 


Here C,, is of the order of 1, w, is close to nz and 6,, is close to zero. 


2 2 
The eigenvalues may be approximated by: k, = K,+ ( zs pe [19] 
t, T t 


In order to derive an improved approximation, we make the assumption: 
K(z)= Ky + €.K,(z) [20] E,(z)= Ey + €-E,(z) [21] 
Here E, and K, are the average values given by [16] and [17], the parameter ¢ is a small 
number and the functions K,(z) and E,(z) are of the order of unity. 
Then the vertical sedimentation-diffusion operator, L, , can be split thus: 
L,=Ty tel; [22] 
Then we split k* = lee + ek [23] and Z(z) =Z, (z)+ é.Z, (z) [24], where Z, and ky are 
given by [18] and [19] if sedimentation is slow enough. 
Equation [15] can be approximated, up to the first order iné , by these two equations: 
AlZo] =O [25] AZ] = g,(z) [26] 
By definition: A[Z,.]=—-Z |Z ]+ e.KoZo—ko PZ) =9 AlZ,]=-£,[Z,]+ e.Ky.Z,—ko’ PZ, 
gi(z)= L, [Zo] — p.K, Zo + pk," Zo 
The solution of [25] is known, as well ask,~. From the theorem of alternatives (Bender 
and Orzag, 1999) we obtain an equation for k,”, so that g,(z) is fully determined. 
Then we can solve [26] in order to derive an expression for Z, . 
Substituting Z,,(z)=Zo,,(z)+¢Z,,,(z) in [14] we can get a better analytical estimation 


ofk,’. 


n 
(D) Discussion and conclusions 


From [12],the quotient of the concentrations of bacteria c, and tracer c; does not 
depend of the horizontal position (x,y), whenever both concentrations are taken in the 
same instant and in the same point, because y,(r;x, y) is the same for both and disappears. 
When the advection-dispersion process has advanced enough, the populations of bacteria 
and tracer are highly diluted. Then the quantities obtained taking local samples at a 
certain depth and at a certain horizontal position may be too scarce to enable good 
measurement results. So, let us suppose that we take samples in vertical columns, from 


z=H (the surface of free water) to z=H-—h (beinghthe length of a vertical sample 
taken with a suitable sampling tube). 


H H 
From [12] it follows that: fel; x, y, Z)dz = Mo. fels z)dz).Wo(t;x, y) [26] 
H-h H-h 
H Beh H 
From [13] we obtain: feGzlde= Ye f,. [Z,(zhdz [28] 
H-h n=0 H-h 


Forming the quotient of the integrals of the concentrations of bacteria c, and tracer c; 


taken along the vertical coordinate between z=H-—h and z=H we obtain this new and 
better relation to apply in practice: 


H A 00 
Joy (sx, y,z}dz M p- Jeol. zhdz yet Ay 
H-h <h Mo» 


7 a a 7. cae [29] 
fer (t3.x, y, z)dz Mor. fer(. z)dz OF Seta 7 ae 
H-h ay ca 


H H 
Here Ay, =fon- [Zyn(z)dz and Ar, = fry». [Zr,(zjdz are parameters that differ between 
H-h H-h 


bacteria and tracer, because not only the eigenvalues, but also the eigenfunctions may be 
different when one population (bacteria) suffers advection, dispersion, perhaps 
sedimentation, and inactivation, meanwhile the other population (conservative tracer) 
suffers advection and dispersion only. 

From [29] and taking into account that for a conservative tracer k*ro=0 , 
asymptotically when t —>oo it follows that: 


H 
fe, (t3.x, y, z)dz 


log,| 4+ ———_|~-K01 +log,(—™*. [30] 
. Mor Apr 
for (t; x,y, z)dz 
H-h 
ars ere 
In this relation, if we neglect sedimentation, k*».0 = Ko = Z [ K(z)dz is inversely proportional 
0 


to the parameter 7,), measured in field experiments with bacteria and tracers. It is an average 
over depth of the field of local kinetic parameters K(z). 


H 
If sedimentation is significant: k740 = =| K(z)dz+ at = 4 
T 
0 


Ss c 


In absence of transport processes the local concentration would decay according to an 
exponential law after a long enough time interval. 

In general, taking into account possible variations from one point to another one, there is a 
whole field of decaying constants, a field of local values of 7) that would determine the 
local capacity for self-purification of the estuarine waters if the combined advection- 
dispersion process were negligible. 

But as this combined process is not negligible, the description of the capacity for self- 
purification must be done by means of the double set of parametersk,,,, kp, » Ap, and Ar, 


for n=0,1,2,3,... A denumerable set of parameters describe the self-purification process. 

In practice a small number may be used to represent the experimental data after the initial 
stages of the process: k750,k751,k°7.1,A,0/Aro> Api/Aro and Ay, /Apo, applying well 
known methods for parameter identification (Thyn and Zitny,2000) . 


From previous work and experimental results, the following half empirical formula for a 
coarse (and generally conservative) estimation of 7,, in hours was proposed: 


T,-T, 1 Ss 
1 a S 
= eK e 2% 4 10.65)C 1 
a S65 

T,: water temperature in °C, s: salinity in Ya oh: angle of the sun’s direction to the 
horizontal plane, in radians; C: fraction of sky covered by clouds; S': solids in suspension, in 
MSA 5 Ko 20.023h", T,; = 20°C, T,5 =52°C, 85 =120% Sy = 800 ae 
It can be used in a gross estimation of the duration of the field experiments in a given scenario 
or as an input to take into account the bacterial die-off in numerical computations using 
hydrodynamic models. 


The following half empirical formula for K, (z) in the simplified analytical model of the 
TT s(z) 


: Ky (:) Ta eo EK fee) 


stratified water body may be of some interest: 
Too(z) 2.3 


=Kge 


B 


being c= , # an average extinction constant in the water body, s(z) the 
1-01.33). cos? a 

vertical salinity profile and @ is the angle between the Sun direction and the horizontal 

plane as defined before. 


H 
} K,(z).dz 
Neglecting sedimentation: ko. = Ky, = 2 7 = = . From the formula for 7), (z) 
TT, s(z) i 
we obtain: 1A ~Kye 1 eae 80 dz)+ Ki to (-e~*) 
90 2.3 H 9 c 


Appendix: A simple and naif approach to bacteria die-off throughT,, 


If we could follow the trajectory of a mass of bacteria c, (t) ee (0) e **" (Phelps (1944) 


D(t) 


C1l0) 1% 
D(t) 
Concentration of bacteria population c,(t); D(t) dilution factor at time t; K, die-off 


; “2 
parameter; 7, = a 


b 
A conservative tracer of concentration c, (t) following the same trajectory would allow to 


Cy (0) 
c,(0) 


formula for rivers) or equivalently c, (t) = 


cue (0) Then log, ont) 


measure dilution: c,(t)= D(t) axe x—K,t+log, would allow us 


to determine the die-off constant K, . 


In a homogeneous body of water, homogeneous both from the advection-dispersion and 
from the bacterial inactivation point of view, this simple approach can be applied. 
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